Absence of hole pairing in a simple t-J model on the Shastry-Sutherland lattice 
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The Shastry-Sutherland model is a two-dimensional frustrated spin model whose ground state 
is a spin gap state. We study this model doped with one and two holes on a 32-site lattice using 
exact diagonalization. When t' > 0, we find that the diagonal dimer order that exists at half-filling 
are retained at these moderate doping levels. No other order is found to be favored on doping. 
The holes are strongly repulsive unless the hopping terms are unrealistically small. Therefore, the 
existence of a spin gap at half-filling does not guarantee hole-pairing in the present case. 



The Shastry-Sutherland (SS) model is an exceptional 
example of a two-dimensional frustrated spin system 
that has an exact solution^ Remarkably, it is also an 
excellent theoretical model for the spin gap material 
SrCu2(B03)2>^ The SS model is a Heisenberg model with 
exchange interaction J on a two-dimensional square lat- 
tice frustrated by diagonal bonds J' as shown in Fig. ^ 
When J/ J' is smaller than a critical value of about 
0.677, the ground state is a direct product of orthog- 
onal singlet dimers residing on the J' bonds. This is a 
spin gap state because the lowest spin excitation involves 
turning a singlet dimer into a triplet. Different experi- 
ments have suggested a range of J/ J' for SrCu2(B03)2, 
out of which 0.635 seems to be the optimal one. A great 
deal of work has been devoted to studying the excited 
states of this model as an effort to explain various exper- 
imental results. It is found that the almost localized na- 
ture of spin excitations leads to superstructures that give 
rise to plateaux in the magnetization curve. Besides its 
importance as a theoretical model for SrCu2(B03)2, the 
SS model possesses a quantum phase transition. One be- 
lieves that there is at least one intermediate state between 
the diagonal dimer state at small J and the Neel state 
at large J. While the nature of the intermediate state is 
still controversial, it seems like the plaquette resonating- 
valence-bond (RVB) state and the helical state are among 
the most likely candidates. Although discussions on the 
quantum phase transition are theoretical, they are not 
totally irrelevant to SrCu2(B03)2. Since its J/ J' is not 
too far from the critical value, it may be possible to shift 
it even closer by applying pressure, substitution, etc. to 

SrCu2(B03)2. 

The relation between disordered spin liquids with a gap 
in the spin excitation and superconductivity has aroused 
a lot of interest. It has been suggested that doping a spin 
gap system may lead to hole-pairing and superconductiv- 
ity. It is therefore interesting and important to identify 
spin gap materials that are Mott insulators with no long- 
range spin order. SrCu2(B03)2 is one such compound. 
Although this compound has not been doped, there is 
no shortage of suggestions that on hopping the SS model 
may exhibit superconductivity. Different types of super- 
conducting states have been suggested at different doping 
levels.-^'*^ One can draw an analogy to a doped two-leg 
spin ladder.^ At half-filling the ground state is a spin 
gap state with singlets forming on the rungs. When two 
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FIG. 1: (a) A 32-site Shastry-Sutherland lattice with peri- 
odic boundary conditions. Also shown are the two types of 
spin interactions and hole hopping terms. I± are reflections 
about the two axis as shown, and [/ is a four-fold rotation 
about the point indicated by •. (b) The corresponding dis- 
tinct reciprocal space vectors. 



holes are introduced, they tend to stay on the same rung 
in order to avoid breaking more singlets. This is a beau- 
tiful example of a spin gap leading to a mechanism for 
hole pairing. The diagonal dimer state of the SS model 
is also a spin gap state. But it differs from the two-leg 
spin ladder in that it is a two-dimensional model. It is 
thus interesting to see whether the same spin gap induced 
pairing mechanism is in effect in the SS model. Further- 
more, doping may introduce more interesting physics to 
the spin background. When J/ J' is close to but below the 
critical value, on doping the system may select other or- 
dered (or disordered) states that can better support hole 
motionn^ Since J/ J' for SrCu2(B03)2 is close to the crit- 
ical value, theoretical discussions on the quantum phase 
transition may become even more relevant if it can be 
doped. Here we report our results on the SS model with 
up to two holes on a 32-sitc lattice using numerical diag- 
onalization. 

We start from the t-J model Hamiltonian, 



(1) 

where tij — t and Jy = J when {ij) is a nearest neighbor 
pair, and tij = t' and = J' when {ij) is a diagonal 
pair as shown in Fig.^a), and zero otherwise. Following 
Ref. 1^, we define a = t' jt. Using the large-C/ expansion 
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TABLE L Ground state energies and symmetries of the model 
with different number of holes A^^ at a = 1.25 and J jt — 
0.3. Nb is the number of basis in that particular symmetry 
subspace. 





E/t 


k 1+ 


/- U la 


Nb 





-12.3 


(0,0) + 


+ + + 


4,708,641 


1 


-15.780 034 


(0,0) + 


+ - 


150,297,603 


2 


-19.291015 


(0,0) + 


+ + - 


601,144,932 



of the Hubbard model as a reference, namely J = At'^/U, 
we take J'/ J — o? . Consequently there are only two 
parameters in our model: a and J jt. Since we want to 
study hole motion in a dimerized spin background and to 
mimic the material SrCu2(B03)2, we fix a = 1.25, i.e., 
J I J' 0.64. The ground state at half-filhng is the SS 
dimerized state with energy —{J+J')N/2 where TV is the 
number of sites. This leaves us with only one adjustable 
parameter J /t when holes are added. Without guidance 
from experimental data, we choose J/t = 0.3 which is 
the appropriate value for high Tc materials. It is also 
the same choice as in Ref. 0. Using exact diagonalization 
we solve this model on a 32-site SS lattice with periodic 
boundary conditions as shown in Fig.^a). The size of a 
unit cell is defined to be 1, i.e., nearest neighbor sites are 
at distance 1/2 apart. Since there are four sites in a unit 
cell, an TV-site SS lattice has N/A translation operations. 
The allowed momenta are shown in Fig. ^b). Spin re- 
flection symmetry Jo- can be used to reduce the size of 
the Hilbert space by a factor of 2 when the number of 
spins is even. In addition, the SS lattice possesses C41, 
(or 4mm) point group symmetry.^ Therefore eigenstates 
can have s, p, or d symmetries, just like the case on a 
square lattice. In our implementation, we choose to clas- 
sify eigenstates according to their symmetry properties 
with respect to the two reflection s 1+ and the four-fold 
rotation U as shown in Fig. ^a).[l^ Table HI shows the 
ground state properties of the model with zero, one and 
two holes. 

The one hole ground state has c?-like symmetry as 
shown in Table ^ Fig. |21 shows the spin structure of 
this ground state. It is obvious that the introduction of a 
mobile hole does not alter the dimerized spin background 
much. Most of the diagonal bonds remain as strong sin- 
glets with (Si • Sj) not much weaker than —3/4. The 
most noticeable feature is an almost free five-spin chain 
formed next to the hole. When a spin is removed from 
one dimer, the resulting free spin tends to minimize en- 
ergy by forming a five-spin chain. In this way the system 
avoids the high energy cost of having a free spin and 
at the same time minimizes disruption to the spin back- 
ground. Note that (Si • Sj) along this "chain" are not 
too far from the corresponding values in an isolated five- 
spin chain with the same J/ J', which are —0.6588 and 
—0.2855 respectively. Also note that it can easily lower 
its energy through hopping — the apex of the chain can 
hop to a diagonal site without disrupting the spin back- 




FIG. 2: Spin structure of the one-hole ground state. The 
hole is at the open lattice point at the center of the lattice. 
Thickness of the lines are proportional to the magnitudes of 
(Si ■ Sj), whose values are also shown for strong singlet bonds. 
Solid and broken lines represent negative and positive num- 
bers respectively. Shaded lines correspond to correlations too 
small to be shown in this scale. Their thickness do not repre- 
sent their magnitudes, but are draw to show their signs only. 

ground. The average kinetic energy per bond on the di- 
agonal and nearest neighbor bonds are — 0.06062 i' and 
— 0.04597 i respectively, showing that diagonal hopping 
is more active than nearest neighbor hopping. Further- 
more, this excitation is quite localized. The bandwidth 
is only 0.1611 1, much smaller than the case on a square 
lattice which is about 0.6 i.^^ 

The two-hole ground state wavefunction is s-like as 
shown in Table The two-hole binding energy Eb = 
E2h - 2Eih + Eo is -0.0309 1 at J/t = 0.3. Although Eb 
is negative, its small value makes it doubtful whether it 
really represents hole-binding. Fig. Ola) shows the hole 
correlation function Chh{f) = ((1 — n,.)(l — no)), which is 
the probability of occurrence of that particular hole con- 
figuration. It clearly shows that the holes prefer to stay 
away from each other, i.e., they are repulsive. The con- 
figuration corresponding to hole-pairing due to spin gap 
effect, where the holes reside on the same diagonal bond, 
is indicated by "A" in Fig.l^Ja). Its hole correlation is the 
smallest. The most favorable hole configuration is when 
they are at maximum distance apart as indicated by "B" 
inFig. |3Ia). Fig. ^ shows its spin structure. A striking 
feature is that next to each hole there is an almost free 
five-spin chain which we have seen in the one- hole ground 
state. Fig. |21 The spin correlation Css(r) = (S^ • So) 
shown in Fig. OJb) is also consistent with the existence 
of such spin chain. In the diagonal dimer state at half- 
filling, Css(r) is zero everywhere except when the holes 
are in the same diagonal bond. Cgsif) at this point re- 
mains the largest in Fig. |3Ib), showing that dimer order 
persists in the two-hole ground state. In addition there 
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FIG. 3: (a) Hole correlation Chhir) and (b) spin correlation 
Casir) of the two- hole model. Dotted lines join inequivalent 
points with the same r. In (b) , open and filled circles represent 
positive and negative correlations respectively. 



FIG. 5: (a) Root-mean-square separation of holes in the two- 
hole ground state at a = 1.25 and 2. The broken line indicates 
the value for a pair of uncorrelated holes on this 32-site lattice, 
(b) Kinetic and spin energies at a = 1.25. 




FIG. 4: Same as Fig. |21 except for the two-hole model pro- 
jected to the subspace where the holes are at the two open 
lattice points. 



are four Css {r) at r — 1/2, l/V^, a/5/2, and V2 that have 
appreciable values. They are due to correlations between 
spins in the five-spin chains. For illustration purpose we 
label these four points in Fig.|2Ib) by (a, b) where a and 
b are sites on a five spin chain such as those marked in 
Fig. 21 For example, the positive correlation (2,4) in 
Fig. El^b) is due to spins 2 and 4 in Fig. 0] Once we 
establish the existence of five-spin chains, we can easily 
understand why holes repel. If the holes are well sepa- 
rated such that a five-spin chain can be formed next to 
each hole, then the local spin energy is minimized by get- 
ting rid of the otherwise free spins, and kinetic energy is 
favored by the hopping motion of the apex of the spin 
chains. 

At t = the two-hole ground state is the exact dimer- 
ized state just like at half-filling but with one dimer re- 



moved. Therefore hole pairing must exist in this limit. 
Let us look at how it is destroyed as t increases. In 
Fig. inja) we plot the root-mean-square separation of 
holes Trms = a/ (r^) at different t/J. The horizontal bro- 
ken line shows the value for two uncorrelated holes on 
this 32-site lattice, which is 1.1914. Let us consider the 
result at a = 1.25 first. At t/J = 0, the holes are at 
the same diagonal bond and r^-^s — 1/V2. But as t/J 
increases, rj-ms increases very fast and at i/J ^ 0.1 it al- 
ready exceeds the uncorrelated value. When t/J > 0.2, 
Ti-ms becomes saturated and the holes become strongly 
repulsive. In Fig. |3Jb) we plot the kinetic and spin en- 
ergies which are the ground state expectation of the first 
and second terms of H in Eq.lQJ respectively. Except 
at very small t/J where r^ms is still increasing rapidly, 
the kinetic energy decreases linearly with t/J while the 
spin energy barely increases. This shows that once the 
system enters the hole-repelling regime, kinetic energy 
of the holes dominates and their motion do not cause 
much disruption to the spin background. This is in total 
agreement with our picture discussed above. 

For the purpose of comparison we also show the cor- 
responding result at a = 2 in Fig. Ela). Note that this 
corresponds to J/ J' — 0.25, which is far from the realis- 
tic parameter for SrCu2 (603)2- With a much larger spin 
gap — 3.7172J compared to 0.6015J at a = 1.25, the 
resulting hole-binding effect persists to larger t/J, and 
Trms reaches the uncorrelated value when t/J is slightly 
larger than 1. However, this represents only a qualitative 
change. The kinetic effect eventually wins and the holes 
are unbound when t/J > 1 . 

So far our results show that diagonal dimer order is 
retained on doping. Next we try to look for "hidden" 
orders that may be associated with charge motion. Our 
previous studj*i^ has shown that the d-wave state of the 
two-hole t-J model on a square lattice has a staggered 
pattern in the current correlation {jkijmn), where the 
current on a bond linking sites k and / is defined by jki = 
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FIG. 6: The current correlation {jkijmn) in the two-hole 
model. The reference bond is denoted by a shaded line with an 
arrow near the center of the diagram. Thickness of solid lines 
are proportional to the magnitudes of the current correlations. 
Thin shaded lines represent correlations that are too small to 
be shown in this scale. 

^tkiicl^cia — cj^Cka)- Tliis Staggered current pattern is 
closely related to the staggered-flux phase based on the 
SU(2) theory. On the SS lattice, a similar staggered- 
flux phase has been suggested by a mean-field studji^. 
However, we know that on a square lattice the staggered 
current pattern exists only in the d-wave state where the 
holes are attractive. But in the present case the ground 
state has s, not d symmetry, and the holes are repulsive. 
Hence we do not expect to see a staggered current pattern 
in this ground state. Fig. clearly shows no staggered 
pattern in the current correlation. Therefore we do not 
find evidence for a staggered-fiux phase in this model. 

Not only does dimerized order persists on doping, we 
believe that doping in fact stabilize it near the quan- 
tum phase transition point. To illustrate this we re- 
peat our calculations at a = 1.2, i.e., J/ J' = 0.694 
which is beyond the quantum phase transition point. The 
ground state at half-filling on a 32-site SS lattice shows 
plaquette order in the dimer correlation^ However, on 
doping with one and two holes the system dramatically 
change back to diagonal dimer order. The spin struc- 
tures and correlations Chhif) and Css(r) are qualitatively 
the same as in Figs. 12 E] and 13 respectively. The one- 
and two- hole ground state energies are — 15. 2 14 070 1 and 
-18.736 734 1, giving a binding energy Eb = -0.0206 1, 
compared to — 0.0309 i at a = 1.25. rmis is 1.2822, com- 
pared to 1.2800 at a = 1.25. Since all results are so sim- 
ilar, our previous conclusions apply to the present case 



with a = 1.2, i.e., the holes are repulsive and the spin 
structure shows strong diagonal dimer order. Since the 
plaquette RVB state exists at most in a very small range 
of J/ J', the fact that on doping it change back to diag- 
onal dimer order means that doping destroys plaquette 
RVB order and stabilize diagonal dimer order. 

To summarize, wc find that doping the dimerized 
ground state of the SS model does not favor any other 
ordered or disordered state. In fact doping favors diag- 
onal dimer order. When a mobile hole is introduced to 
the dimerized state, the "free spin" that would have been 
formed minimize energy by forming a five-spin chain with 
two of its nearest dimers. This excitation is localized and 
has a small dispersion compared to the t-J model on a 
square lattice. When more than one hole is present, their 
respective spin chains tend to avoid each other. This give 
rise to short-range repulsion between the holes. Although 
there may exist long-range attraction between the holes 
mediated through the chains. Fig. |3Jb) shows that the 
spin energy is not sensitive to changing spin configura- 
tions. Therefore we consider long-range attraction un- 
likely. There is another competing short-range pairing 
mechanism — the spin gap naturally leads to pairing of 
holes on the same diagonal bond when t — 0. But we find 
that this pairing mechanism is in effect only when t/ J is 
smaller than about 0.2 at a = 1.25. Note that this value 
corresponds to U/t ~ 4t/J = 0.8, which is far from the 
strongly interacting regime of the Hubbard model. As a 
result, the idea of treating the dimerized ground state of 
SrCu2(B03)2 as a Mott insulating state of an underly- 
ing Hubbard model will fail. Consequently we consider 
t/J < 0.2 too small to be realistic. The implication of 
our result is that the existence of spin gap in a disordered 
spin liquid state does not guarantee charge binding. 

Finally, we remark that our results do not completely 
rule out hole pairing in the SS model because we have 
studied the case t' > only. Just like in the case of a 
square lattice, the sign of t is irrelevant. On a square 
lattice, t' breaks electron-hole symmetry and Neel order 
of the spin background is better preserved with t' > 0. 
This property is consistent with the phase diagram of 
electron-doped high Tc superconductors>i^ Similarly, we 
have shown that on a SS lattice diagonal dimer order of 
the spin background persists on doping when t' > 0. But 
without experimental results we cannot decide whether 
t' > corresponds to electron-doping in the present case. 
Nevertheless, it is meaningful to study the case t' < 
where diagonal hopping is discouraged. These works are 
in progress and will be presented elsewhere. 
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